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Abstract. It has long been known that a large fraction of disc galaxies are lopsided. We simulate 
three different mechanisms that can induce lopsidedness: flyby interactions, gas accretion from 
cosmological filaments and ram pressure from the intergalactic medium. Comparing the mor- 
phologies, HI spectrum and m = 1 Fourier components, we find that all of these mechanisms can 
induce lopsidedness in the gaseous component of disc galaxies. In particular, we estimate that 
flybys can contribute to ~ 20 per cent of lopsided galaxies. We compare our simulations with 
the observations of NGC 891, a lopsided, edge-on galaxy with a nearby companion (UGC 1807). 
We find that the main properties of NGC 891 favour a flyby event for the origin of lopsidedness 
in this galaxy. 
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1. Introduction 

A high fraction of disc galaxies are lopsided, i.e. their gas and/or stellar component 
extend further out on one side of the galaxy than on the other (Baldwin, Lynden-Bell & 
Sancisi 1980; Richter & Sancisi 1994; Rix & Zaritsky 1995; Zaritsky & Rix 1997; Swaters 
et al. 1999; Bournaud et al. 2005; see Sancisi et al. 2008 for a review). 

Richter & Sancisi (1994) show that the lopsidedness of a galaxy can be inferred from 
asymmetries in its global HI profile, and estimate, from the analysis of 1700 HI spectra, 
that > 50 per cent of disc galaxies are lopsided in the gaseous component. The kinematics 
of the gas is often affected by lopsidedness: Swaters et al. (1999) find that the rotation 
curve of lopsided galaxies is rising more steeply on one side than on the other. Rix & 
Zaritsky (1995) and Zaritsky & Rix 1997, using near-infrared photometry of nearly face- 
on spiral galaxies, show that even the stellar component is lopsided in ^ 30 per cent of 
their sample. 

The hypothesis that lopsidedness is due to galaxy interactions has been long discussed. 
Based on optical images, Odewahn (1994) finds that 71 of 75 lopsided Magellanic spirals 
have a nearby companion. However, Wilcots & Prescott (2004) obtain HI data of 13 
galaxies from Odewahn (1994) and show that only four of them have Hl-detected neigh- 
bours. Thus, either lopsidedness is not related to galaxy interactions, or the asymmetries 
produced by these interactions are long-lived (surviving for > 2 orbital times after the 
encounter) and the lopsidedness persists even when the companion is quite far-off. 

From the theoretical point of view, the N-body simulations by Walker, Mihos & Hern- 
quist (1996) suggest that minor mergers can induce lopsidedness over a long timescale 
(> 1 Gyr). However, Bournaud et al. (2005) indicate that the lopsidedness produced by 
minor mergers disappears when the companion is completely disrupted. Since most of 
observed lopsided galaxies are not undergoing mergers, the minor- merger scenario does 
not seem viable. Bournaud et al. (2005) indicate that the most likely mechanism to pro- 
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duce lopsidedness is the accretion of gas from cosmological filaments. Alternative models 
suggest that baryonic lopsidedness can be induced by a lopsided dark matter halo (Jog 
1997) or by the fact that the disc is off-centre with respect to the dark matter halo 
(Levine & Sparke 1998). 

Here we address the problem of the origin of lopsidedness by means of A^— body/smoothed 
particle hydrodynamics (SPH) simulations. In particular, we re-analyze in more detail 
the hypothesis of gas accretion, already proposed by Bournaud et al. (2005), and we con- 
sider two new possible scenarios: the role of flyby interactions with smaller companions 
and that of ram pressure from the intergalactic medium (IGM). 

2. Methods 

We consider three different processes: flybys (i.e. interactions between a target and an 
intruder galaxy which do not imply merger), gas accretion from cosmological filaments 
(Keres et al. 2005) and ram pressure from low-density IGM. For each of these cases we 
make a large number of simulations with good spatial and mass resolution (the softening 




Figure 1. Projected density of gas in the simulations. Top (Bottom) panels: the target galaxy 
is seen edge-on (face-on). From left to right: flyby, gas accretion and ram pressure scenario. 
The left-hand panels measure 80 kpc per edge. The central and the right-hand panel measure 
100 kpc per edge. The time elapsed from the beginning of the simulation is 0.64, 0.72 and 
1.0 Gyr in the flyby, gas accretion and ram pressure case, respectively. The density goes from 
2.23 X 10"^ to 2.23 x 10^ Mq pc"^ from 5.60 x 10"^ to 2.23 x 10^ M© pc"^ and from 7.05 x 10"^ 
to 2.23 X 10^ Mq pc~^, in logarithmic scale, in the flyby, gas accretion and ram pressure case, 
respectively. 
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length and the mass of each particle are 0.1 kpc and 2.23 x 10'* Mq for gaseous and stellar 
particles, and 0.2 kpc and 2.23 x 10^ Mq for dark matter particles, respectively). The 
A^-body/SPH code used is GASOLINE (Wadsley, Stadel & Quinn 2004). The details of 
the simulations have been described in Mapelli et al. (2008a). Here we provide only the 
fundamental characteristics. 

Flybys. We make different runs, with different mass ratios between the intruder and 
the target (1:20 to 1:2), with different masses of the intruder and of the target, with 
different impact parameter and relative velocity. 

Gas accretion. The simulated galaxy accretes a cylinder of uniform gas with relative 
velocity of 100 km s~* and with mass accretion M = 2 — 6Mq. Different runs have 
been made with different inclinations of the gas cylinder with respect to the plane of 
the galaxy. This model, although very simplified, is supported by the gaseous filaments 
observed in cosmological simulations (Keres et al. 2005). 

Ram pressure. The determination of the density of IGM in the field and in poor groups 
is currently an open issue. Recent observations (Pisano et al. 2004; Freeland, Cardoso & 
Wilcots 2008) suggest that the IGM density in groups is > 10~^ cm~'^. At such densities, 
relative velocities > 100 km s^^ should lead to ram-pressure stripping (Gunn & Gott 
1972). Thus, we simulate a galaxy moving with a velocity of 200 km in a uniform 
IGM with density 5 x 10"^ cm'^. 

3. Results 

The main result of our simulations is that all these mechanisms (flybys, gas accre- 
tion and ram pressure) produce lopsidedness in the gaseous component of the simulated 
galaxy. This can be seen qualitatively by the density maps shown in Figure 1. Such re- 
sult is important, as it indicates that lopsidedness is not due to a single process, but 
may originate from different mechanisms. The total stellar component of the galaxy is 
significantly lopsided only in the case of flybys (upper panels of Figure 2). Interestingly, 




Figure 2. Projected density of stars in the simulations. Top panels: all stars. Bottom pan- 
els: young stars (< 100 Myr). From left to right: flyby, gas accretion and ram pressure sce- 
nario. The left-hand panels measure 80 kpc along the long side. The central and the right- 
hand panel measure 100 kpc along the long side. The time elapsed from the beginning of the 
simulation is 0.64, 0.72 and 1.0 Gyr in the flyby, gas accretion and ram pressure case, respec- 
tively. Top (bottom) panels: the density goes from 2.23 x 10"^ to 2.23 x 10^ Mq pc"^ (from 
7.05 X 10"® to 2.23 x 10^ M@ pc"^), from 7.05 x 10"^ to 2.23 x 10^ Mq pc"^ (from 7.05 x lO"** 
to 2.23 X 10^ Mq pc"^) and from 5.60 x 10"^ to 2.23 x 10^ Mq pc"^ (from 7.05 x 10"'^ to 
2.23 X 10^ Mq pc~^), in logarithmic scale, in the flyby, gas accretion and ram pressure case, 
respectively. 
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the young stellar component (i.e. the stellar particles with an age < 100 Myr) is strongly 
lopsided also in the gas accretion scenario (bottom panel of Figure 2). This indicates 
that, in the gas accretion scenario, the lopsidedness of the stars may be induced by star 
formation in a lopsided gaseous component. A more quantitative estimate of lopsidedness 
is given by the normalized strength of the Fourier component m = 1 (Al) as a function 
of radius, shown in the left-hand panel of Figure 3. The Fourier analysis confirms that all 
the three mechanisms induce lopsidedness in the gaseous component. The flyby scenario 
produces lopsidedness both in the old and in the young stars, the gas accretion scenario 
only in the young stars, whereas the ram-pressure model does not induce lopsidedness in 
the stars. 

We also analyze the kinematic properties (HI spectrum and rotation curve) of the 
simulated galaxies. The HI spectrum is shown in the right-hand panel of Figure 3. The 
simulations (blue line) are compared with the data of the lopsided galaxy NGC 891 (green 
circles, from Richter & Sancisi 1994). The main feature of lopsidedness in the HI spectra 
is the asymmetry of the two peaks of the spectrum. Our simulated galaxies show only a 
moderate lopsidedness in the HI spectra. In particular, the flyby and the ram-pressure 
case agree very well with the data of NGC 891. 

3.1. Flyby rate 

We have shown that flybys can induce lopsidedness in galaxies, both in the gaseous and 
in the stellar population. We now derive which fraction of disc galaxies may be lopsided 
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Figure 3. Left-hand panel: Normalized strength of the Fourier component m = 1 as a function 
of radius, for gas (top panel), global stellar content (central panel) and young stellar popula- 
tion (bottom panel). Open circles connected by red solid line: flyby scenario; open triangles 
connected by blue dotted line: gas accretion; open squares connected by green dashed line: ram 
pressure. Right-hand panel: HI spectra of the simulated galaxies compared with the observations 
of NGC 891. Open green circles: data from Richter & Sancisi (1994). Solid blue line: simulations. 
Top panel: flyby interaction; central panel: gas accretion; bottom panel: ram pressure. 
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due to flybys (an analogous estimate is currently hard to do for gas accretion and ram 
pressure, given the insufficient observational constraints). The density of galaxies which 
are lopsided as a consequence of flybys may be derived from the density of coUisional 
ring galaxies (CRGs, see Mapelli et al. 2008b and references therein) in the local Uni- 
verse (ncRG — 5.4 X 10~^h^Mpc~^, where h is the Hubble parameter, Few & Madore 
1986). In fact, the galaxy interactions which induce lopsidedness are similar to the ones 
which produce the CRGs, although with some crucial difference (e.g. the range of impact 
parameters, the inclination angles, the lifetimes). Taking into account this analogy, the 
density of galaxies which are lopsided due to a fly by (nfly) can be derived as: 
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where 5max ^ '^Rd and &CRG ^ 0.36 i?^ are the maximum impact parameter to induce 
lopsidedness and to form a CRG, respectively {Rd is the disc scale- length), whereas 
t\o-p and tcKG are the lifetime of lopsidedness and of a CRG, respectively. The term 
(1 — cos 45°) accounts for the fact that only interactions with inclination angle ^ 45° 
produce CRGs (Few & Madore 1986). Assuming t\op ^ 1 Gyr and ^crg ^ 0-5 Gyr, 
as suggested by simulations, we obtain nfly ^ 4.5 x 10~^ Mpc~^. Since the density 
of relatively bright disc galaxies (i.e. with absolute magnitude ^ — 19 mag) is ridisc ~ 
4.9 X 10~^ Mpc^'^ (Few & Madore 1986), this suggests that only one galaxy every 
~ 11 spiral galaxies can be lopsided as a consequence of flybys. As lopsidedness of the 
gas component is observed in ~ 50 per cent of disc galaxies (Richter & Sancisi 1994), 
flybys can reasonably contribute to ~ 20 per cent of lopsided galaxies. 

3.2. Comparison with NGC 891 

Finally, we compare the properties of the simulated galaxy with those of one of the most 
studied lopsided galaxies, NGC 891 (Oosterloo, Fraternali & Sancisi 2007 and references 
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Figure 4. HI spectrum of the simulated intruder galaxy (i.e. the small companion) in the 

flyby scenario. 
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therein). This galaxy is mildly lopsided (right-hand panel of Figure 3) and has some 
interesting features. In particular, NGC 891 has a conspicuous HI halo 30 per cent of 
the total observed HI). Furthermore, the observations show a gaseous filament extending 
up to ~ 20 kpc vertically from the disc and located at ~ 10 kpc from the centre of the 
galaxy. Finally, NGC 891 has a smaller, gas-rich companion, UGC 1807, located at a 
projected distance of ~ 80 kpc, in the direction of the above mentioned gaseous filament. 

The simulations reproduce quite well the lopsidedness of NGC 89l[|. In particular, the 
flyby and the ram pressure scenarios match the HI spectrum (right-hand panel of Figure 
3). However, the gas accretion and the ram pressure simulations cannot reproduce the 
observed filament. Only the flyby scenario produces a sort of gaseous filament between 
the target and the intruder (left-hand panel of Figure 1). This filament is similar to the 
one observed in NGC 891, although there are some differences (e.g. the simulated fila- 
ment extends with almost uniform density between the target and the intruder, whereas 
the observed filament is denser in the region which is close to NGC 891). Thus, our 
simulations suggest that the lopsidedness of NGC 891 may be due to fiybys. However, 
further comparisons with data are required. For example, in Figure 4 the HI spectrum of 
the simulated intruder galaxy is shown. This spectrum appears mildly lopsided, similarly 
to that of the target galaxy. It would be interesting to look at the spectrum of UGC 1807 
(and, in general, of the companions of other lopsided galaxies) and to see if they appear 
lopsided, too. 
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